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Abstract: In this paper we address the problem of detection of entanglement 
using only few local measurements when some knowledge about the state is 
given. The idea is based on an optimized decomposition of witness operators 
into local operators. We discuss two possible ways of optimizing this local 
decomposition. We present several analytical results and estimates for optimized 
Qh. detection strategies for NPT states of 2 x 2 and N xM systems, entangled states 

in 3 qubit systems, and bound entangled states in 3 x 3 and 2x4 systems. 

& . 

1 Introduction 

One of the key problems in experiments in quantum information theory is the 
k>( , generation and detection of entanglement [TJ. The generation of entangled states 

; i ■ is practically always aimed at particular states that can then be used for var- 

ious applications in quantum information processing. Very often one wants to 
produce certain pure states. The generation process in a laboratory is, however, 
never free of imperfections and noise. The states that one produces may, but 
not necessarily have the desired properties. In particular they may, but do not 
have to be entangled. For this reason it is important to develop efficient and 
easy to apply experimental procedures to detect entanglement. 

Obviously, the ultimate goal of detection is to characterize the entanglement 
quantitatively, and to identify the regions in the parameter space which would 
allow to maximize entanglement useful for the particular application. Before 
this ambitious goal is realized, however, it is important to know whether the 
state one is dealing with is entangled at all, or not. The reasons for that are 
at least threefold: i) it is interesting from the fundamental point of view; ii) 
it is important to know it before applying more complex experimental tools 
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to quantify entanglement; iii) it is essential, if one wants to use the generated 
states for any of the distillation or purification protocols [2] ■ 

For the detection of entanglement several strategies are known: First of all 
there is the possibility of quantum state tomography P] , and then direct appli- 
cation of known necessary or sufficient entanglement criteria 0]. Determining 
the density matrix via quantum state tomography requires, however, typically 
a lot of measurements. Note also that only for 2 x 2 and 2x3 systems nec- 
essary and sufficient criterion is known, the famous Peres-Horodccki criterion 
of positivity of partial transpose j^j. In general, it might be difficult to find a 
sufficient criterion for a given state. One could also look for a violation of some 
Bell inequalities, although there exist entangled states (bound entangled states) 
that do not violate any known Bell- like inequality 0. In fact, there is even a 
conjecture that these states admit a local hidden variable model [7]. Looking 
for a violation of Bell inequalities may thus not be sufficient. 

Recently, there have been several proposals for the detection of entanglement 
without estimating the whole density matrix [HI EH ■ Although being attrac- 
tive from both experimental and theoretical point of view, these possibilities 
have some disadvantages. The recent proposals require collective measurements 
on several qubits or the construction of quantum gates and networks. These 
requirements arc possible, but not easy to fulfill with the present experimental 
techniques of local measurements. 

On the other hand, the methods for the detection of entanglement mentioned 
above are in some sense too general for experiments. They assume that no a 
priori knowledge about the density matrix is given. But, as we already men- 
tioned, in an experimental situation one usually tries to prepare some particular 
state. Although one fights with the problems of noise and imperfections, in any 
case the produced state cannot be considered completely arbitrary. 

In this paper, which expands our earlier study |1 1| . we address the problem 
of checking whether a state g is entangled or not, when some knowledge about 
the density matrix is given. We want to solve this problem using only local 
von Neumann measurements which can be implemented in a laboratory using 
the present techniques. We also aim at using the smallest possible number of 
measurements. 

The scheme we use for the detection of entanglement relies on the well known 
concept of witness operators E| • Let us briefly recall what these operators 
are. A Hermitean operator W is called an entanglement witness detecting the 
entangled state g e if Tr(W g e ) < and Tr(W g s ) > for all separable states 
g s . So, if we have a state g and we measure Tr(Wg) < 0, we can be sure that 
g is entangled. For every entangled state there exists an entanglement witness. 
There has also been an enormous progress in constructing witnesses for different 
classes of states 

After having constructed a witness operator we decompose it into a weighted 
sum of projectors onto product vectors. In this way the expectation value of 
the witness can be measured locally: Alice and Bob measure the expectation 
value of the projector and add their results with the weights to receive the 
expectation value of the witness. Knowing the value of Tr(W g) they can decide 
if the state is detected to be entangled, or not. The decomposition of the witness 
into projectors onto product vectors should be optimal, i.e. contain in a certain 
sense a possibly small number of terms. 

The paper is divided in four sections. In the first section we illustrate our 
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proposal and our notation with a simple example: We consider a pure two-qubit 
state affected by white noise. The main parts of this paper are the extensions of 
this example: In the second section we investigate the two-qubit scenario if the 
noise is not white. We also investigate the limits of our scheme and the errors 
that may occur. The third section deals with extensions of the decomposition 
of the witness from the first section. We derive decompositions for witnesses 
in three-qubit case. With these operators GHZ-states and W-states can be 
detected. Eventually, we discuss how one can decompose witnesses for bipartite 
systems in higher dimensions. In the last section we apply our method to bound 
entangled states. We mainly discuss how some types of bound entanglement in 
3x3 and 2x4 systems can be detected. 

2 Two qubits with white noise 

We consider a setup that is intended to produce a particular pure state |V')(V , Ij 
but due to the imperfections some noise is added. So, it produces a mixed state 
g of the form 

g(p,d) :=p^)(Vl + (l-pV, (1) 

where we know that the noise added to the state is close to the totally mixed 
state, i.e. 

\W-\n<d. (2) 

If d = the noise would be white, but in general this docs not have to be the 
case. In the beginning, we do not make any assumption on but when we 
discuss the case d > we restrict ourselves to the important case that is a 
Bell state: \ip+) = l/\/2(|01) + |10>). Further we do not make any restriction 
on p, i.e. we assume that it is unknown to us. As a byproduct the witness 
will later donate us a possibility to determine p. Our aim is to give a scheme to 
check whether g(p, d) is entangled or not. 

Now we have to clarify our notation. We will denote the set of all density 
matrices by M and the set of separable matrices by S. We can endow the 
space of all matrices with the scalar product (A\B)hs '■= Tr(A(B^)) and the 
corresponding Hilbcrt- Schmidt norm \\A\\ := y/Tr(A(A^)). For any g £ M and 
r > we can define by 

B{g,r):={g' £M,\\g' -g\\<r] (3) 

the ball around g with radius r. We denote by B p ^ the ball which must include 
g(p,d) for a given p and d. It is given by B p ^ = 5(^-^1 + p\ip + )(ip + \, (1 — 
p)d). Finally we define L as the line between 1/4 and we have: L = 

{g(p,0), p£ [0,1]}. 

2.1 Construction of the witness 

The optimal witness for an entangled g(j>, 0) is easy to construct |13| . First, 
one has to compute the eigenvector corresponding to the negative eigenvalue of 
g(p, 0) Tn , then the witness is given by the partially transposed projector onto 
this eigenvector. 

If the Schmidt decomposition of \%p) is \ip) = a\01) + b\l0) with a, b > 0, the 
spectrum of g(p, 0) Tb is given by {(1 —p)jA. + pa 2 ; (1 — p)/4 + pb 2 ; (1 — p)/4 + 
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pab; (1 - p)/4 - pab}. Therefore g(p, 0) is entangled iff p > 1/(1 + Aab). The 
eigenvector corresponding to the minimal eigenvalue A_ is given by 



= -=(|00) - |11», 



(4) 



and thus the witness Wo is given by 
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Note that this witness docs neither depend on p, nor on the Schmidt coefficients 
a,b. It detects g(p,0) iff it is entangled, since we have Tr(|0_)(0_ \ Tb g(p, 0)) = 
Tr(\(j)-)((j)-\g(p, 0) Tb ) = A_. For our special case of white noise Wq has also 
the advantage that from Tr(Wog(p 7 0)) > it follows that g(p, 0) is separable. 
This is not a general property of witnesses, and indeed if the noise is not white 
this is not true anymore. 

2.2 Decomposition of the witness 

For an experimental setup it is necessary to decompose the witness into op- 
erators which can be measured locally. Thus we need a decomposition into 
projectors onto product vectors of the form 



Such a decomposition can be measured locally: Alice and Bob measure the 
expectation value of the |ei) (e^ | <g> \fi)(fi\ and add their results with the weights 
Ci. One can construct such a decomposition in many ways, but it is reasonable to 
do it in a way which corresponds to expenses of Alice and Bob which are as small 
as possible. There are several possibilities to define an optimal decomposition: 

One possibility is to look for the optimal number of product vectors (ONP), 
i.e. one can try to minimize k in ©. This optimization strategy looks very 
natural and has already been considered in the literature. For 2x2 systems it 
was proven in ^1] that in general the ONP, which we denote as is five. Also 
a constructive way for computing this optimal decomposition was given. 

What is the "cost" Alice and Bob have to pay when measuring W via such a 
decomposition? It is the number of measurements they have to perform. When 
we talk about measurements here, we consider only von Neumann measure- 
ments. We do not look at POVMs since their implementation would require 
additional ancilla systems. One measurement on Alice's side in this sense con- 
sists of a choice of one orthonormal basis for Alice's Hilbert space. For an 
particle with spin s one may interpret this as the choice of a direction for a 
Stern-Gerlach-like apparatus. Alice sets up her device in the desired direction 
and is able to distinguish between 2s + 1 different states. 

For one local measurement Bob also has to choose an orthonormal basis 
in his Hilbert space; all together this yields one orthonormal product basis for 




(6) 
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both. Thus, if we are in a TV X N system a term of the form 

N 

]T c kl \A k ){A h \®\B l ){B l \ (7) 

with (A s |j4 f ) = (B s \B t ) = 6 s t can be measured with one collective setting 
of measurement devices of Alice and Bob. Alice and Bob can discriminate 
between the states \A k Bi), measure the probabilities of these states and add 
their results with the weights c k i using one collective setting and some classical 
communication. We call such a collective setting of measurement devices a local 
von Neumann measurement (LvNM). 

It is therefore reasonable to find a decomposition of the form 

m N 

W = Y J Y. c ^)(A k \®\Bl){B\\ (8) 

i=i k,i=i 

with (A\ \A\) = (B\\B\) = S st and an optimal number of devices' settings (ONS), 
i.e. a minimal m. In this sense m is the minimal number of measurements Alice 
and Bob have to perform. The construction of a decomposition of the form JHJ, 
and the determination of the minimal m is the problem we want to solve in this 
section. 

Please note that a decomposition like © with the minimal fc_ (ONP) re- 
quires in general fc_ LvNMs because the out-coming vectors on Alice's side 
|ej)(ej| do not have to be orthogonal. 

We also would like to emphasize that a decomposition of the form (JSJ is 
more general than a decomposition into a sum of tensor products of operators: 

rn 

W = J2-7iAi®Bi. (9) 
i=i 

The decomposition © has the advantage that Alice and Bob do not have to 
distinguish between some states, they only have to measure locally some expec- 
tation values of Hermitean operators. A decomposition like (JSJ can be written 
in the form of if for all i the matrices {c l kl ) are of rank one. In the follow- 
ing we will see that for qubit systems there is not a big difference between © 
and @. From the optimal decomposition in the sense of JHJ we can derive a 
decomposition of the form © where some of the operators are the identity (1), 
so they do not require new measurement settings. For N x A^-systems we will 
see that it is straightforward to derive the optimal decomposition in the sense 
of ©. 

Our witness Wo has a special form: It is a partially transposed projector 
(ip\ TB . It is therefore natural to look at the Schmidt decomposition of |-y>) = 

a|00) + For our Wo we have the special case a = 1/^/2 = —f3, but we 

want to deal with the most general \tp). 

When we compute the ONP-decomposition with the minimal /c_ according 

to we arrive at 

mW T B = (* + P)lj2\A' i B>)(A>B>\ - a/3(|01)(01| + |10)(10|), (10) 
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where we have used the definitions 



\A[) = cos(0)|O) + e - ** sin(0)|l) = \B[) 

\A' 2 ) = e _i * cos(6»)|0) + e<* sin(6»)|l) = 

|A 3 ) = |A' 1 ) + |A' 2 ) = |^) 

cos(#) = y/a/(a + /3) 

sm(0) = y//3/(a + l3). (11) 

This decomposition into five product vectors requires four correlated settings for 
Alice and Bob. But we can measure Wq with less settings: If we define the spin 
directions by \z+) = |0>, \z~) = |1), |x±) = ^(|0) ± |1), |y±) = ^(|0) ±i\l) we 
have the decomposition 

\^)0j\ Tb = a 2 \z + z + )(z + z + \ + 2 \z~z^){z~z-\ + a(3(\x+x + ){x + x + \ + 

+\x~x~)(x~x'\ - \y + y~)(y + y~\ - \y~y + ){y~y + \) 

= - (l ® 1 + a z ®a z + (a 2 - (3 2 )(a z ® 1 + 1 ® cr z ) 

+2a/3(cr x ®cr :!; +cry(8)cr y )). (12) 

This decomposition into six product vectors requires only a measurement of 
three settings: Alice and Bob have only to set up their Stern-Gerlach devices in 
the x-, y- and z-direction to measure \ip)(ip\ TB . 

Now we want to prove that three LvNM are really necessary. Our proof is a 
special case of a theorem about N x N systems we will show later. But in the 
two-qubit case the proof is particularly simple, and therefore we present it here 
separately. 

Proposition 1. In a two qubit system a decomposition of \ip)(ip\ TB of the 

form (JSJ) requires at least three measurements. 

Proof. Consider a decomposition requiring two measurements: 

2 2 

Wm T ° = E c iM$)( A t\ ® \b)){b}\ + £ 4K>M?I ® I^><b?|. (is) 

With the help of a Schmidt decomposition as above we can write \ip)(ip\ TB = 
2, 7 =o °i ® °j witn 
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Note that the 3x3 submatrix in the right bottom corner is of rank 3. Now we 
write any projector on the rhs of 113|) as a vector in the Bloch sphere: \A\) (A\ | = 
£j=o s t a i i s represented by the vector s A i = (1/2, sf-, s 2 , sf) and I^X^I ^y 
s A i = (1/2, — sf , — s 2 , — s^); \B\){B\\ can be written similarly. If we expand 
the first sum on the rhs of (|13|) in the (<Xj <g> <7j) basis, the 3x3 submatrix in the 
right bottom corner is given by (c^ — c\ 2 — c 21 + c\ 2 )(s^ , s; 4 , s 3 4 ) T (sf , s^, S3 3 ). 
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This matrix is of rank one. The corresponding submatrix from the second sum 
on the rhs of H13fl is also of rank one and we arrive at a contradiction: No matrix 
of rank 3 can be written as a sum of two matrices of rank one. □ 
The idea of a generalization of this proposition to N x N systems is straight- 
forward: One expands both sides of an equation of the type l|13|) in some product 
basis of the space of all operators. Then one tries to reach lower bounds for the 
rank of some submatrix of the coefficient matrix l|14|l and upper bounds for the 
rank of the matrix corresponding to one LvNM. This gives a lower bound for 
the number of LvNMs. 



3 The case of non-white noise 

Our investigation of this case proceeds in two steps. First, we argue why the 
witness Wq should also be used in this instance. We also show that under some 
circumstances, i.e. if Tr(Wa@(p, d)) is large enough, we can make a sure decision 
that g(p, d) is separable. For the case that we cannot make a sure decision we 
derive analytical bounds and show numerical estimates for the error. This is 
done in the second step. 

We would like to remind the reader that if d > 0, we only consider maximally 
entangled states, so everywhere in this section we set a = b = 1/V2. 

3.1 Properties of W 

Let us note some typical distances in M: 

Remark 1. (a) The states on dM, the boundary of M, closest to 1/4 have a 
distance of l/\/l2 « 0.29. The states with the largest distance to 1/4 are just 
the pure states, for which the distance is V3/2 w 0.87. 
(b) The states on OS closest to 1/4 have also distance of 1/^12. 
Proof. These distances can be calculated simply by maximizing or minimizing 
|| 1/4 — g\\ under some conditions, for instance for (b) under the condition that 
dct(g Ts ) = 0. □ 
From part (a) it follows that we have to assume that d < 1/ ^/l2. Please note 
also that the point where L crosses the border of separability, dS, is given by 
£>(l/3,0) and this point has just the distance 1/Vl2 from 1/4, this means it is 
as close as possible. 

Remark 2. For all p the point on dS closest to g{p, 0) is given by £>(l/3, 0). 
Proof. For p < 1/3 we have mentioned it already. For p > 1/3 we have a look 
at the set 

Af{W ) := {g G A/; Tr(W g) = 0} . (15) 

We have £>(l/3, 0) £ N(Wq), and one can directly compute that for an arbitrary 
N £ Af{W ) and g(p,0) the relation (#(1/3,0) - g(p, 0)\g(l/3, 0) - N) HS = 
holds. Therefore JV(Wq) and L are orthogonal. From this and the properties of 
a witness operator the claim for p > 1/3 follows. □ 

Remark 3. The witness Wo is the best possible witness in the following sense: 
For all p and d there is no other witness that detects a subset of B Pi d as en- 
tangled which has a bigger volume than the subset of B Pt d which is detected as 
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entangled by Wo- 

Proof. From the second remark it follows that if B p ^ H S = the witness Wo 
detects the whole ball, independent of d. 

Thus, we can assume that p has a value such that B Pi d contains separa- 
ble and entangled states and that there exists a witness W' , which detects a 
bigger volume of B p ^. Then we look at the set Af(W'), defined analogous to 
JTHJ. If p < 1/3 the minimal distance between matrices of N(W') and g(p,0) 
must be smaller than the minimal distance between TV" (Wo) and g(p, 0), which is 
\\g(p, 0) — g(l/3, 0)|| . This means that there exist entangled states (and therefore 
states on dS) which are closer to g(p, 0) than g(l/3, 0) and we have a contradic- 
tion to Remark 2. If p > 1/3 the minimal distance between matrices of Af(W') 
and g(p, 0) must be bigger than the distance between 7V(Wo) and g(p, 0), but 
then W' must "detect" the separable state p(l/3,0), which is a contradiction. 
□ 

The critical reader may ask at this point why we are so innocent and use the 
term of a "volume" in M. It seems that it is difficult to say something about 
volumes in our norm, since we are not in the R™ with the Euclidean norm |15II16| . 
But we are not too far away from the K™ : We can write any density matrix 
as g = 2j=o MiCi where fii G K, Go ~ 1 and the G±, G\5 arc the traceless 
generators of the SU(4). The Gi can be normalized in a way that they form an 
orthonormal basis: (Gi\Gj) hs = $ij ■ We will explain and use this decomposition 
in greater detail later. Here we only point out that \\g\\ 2 = J2i=oA*f an< ^ so our 
norm in the space of Hcrmitcan operators is just the Euclidean norm in R 16 . M 
corresponds to a subset of a 15-dimcnsional hyperplane since /xo is fixed, and 
the usual formulas for volumes of balls can be applied. 

Now having proven that it is reasonable to use Wo, one may ask what the 
expectation value of Wo tells us for the case d > 0. From the definition of a 
witness operator it follows that if the expectation value is negative we can be sure 
that g(p, d) is entangled. On the first view it seems that if it is positive we can 
not make a sure decision if g(p, d) is entangled or not. But if Tr(Wog(p 7 dj) ^> 
this means that g(p, d) is far away from Af(Wo) and since we know that g(p, d) € 
B Pt d for some p, the state g(p, d) must be separable. So there must exist a r(d) 
such that from Tr{g{p 1 d)Wo) > r it follows that g(p,d) is separable. We can 
directly compute r, but first we introduce some new definitions. 

If we measure Tr(gWo) = a the expectation value a tells us that g is in some 
hyperplane cutting M. This hyperplane is orthogonal to L, as we have shown 
in the proof of remark 2. It intersects this line at some g(q, 0) and so we denote 
this plane by V(q). For instance an essential part of the proof of Remark 2 was 
the statement that 7^(1/3) = TV (Wo). The connection between the expectation 
values of Wq and the planes V is given by: 



Since we do not deal with general g, but with some g(p, d), it will be useful to 
consider the set of all possible g(p, d). This is 



Tr(gW ) = a g e V(q) with q = - - -a. 



(16) 




(17) 
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We can now define the intersections 



SP(q) := SnV(q), 
KP(q,d) := K(d)nV(q), 
BP{q lP ,d) := B M nP(g), 

XP(g) := B(^)n%). (18) 

SP is the set of all separable states which belong to one possible expectation 
value. KP is the set of all possible g(p, d) yielding the same expectation value. 
It is clear that BP(q,p,d) C KP(q,d). Since we know already from Remark 1 
that the states in P(l/4, 1/V12) are separable, we can conclude that XP{q) C 
SP(q). For the sake of notational simplicity we often suppress the parameters 
q,p,d. 

The strategy of computing r is now clear: We have to find the values of q for 
which KP(q) C SP(q) holds. For the corresponding expectation values of Wo 
our knowledge that g(p, d) £ K guarantees us that g(p, d) is separable. Since it 
is difficult to characterize SP we replace it by XP. This replacement is justified 
later. 

Proposition 2. Let 

H- i2 -\Z(s- <p )(i-4 (19) 

Then we have: 



Tr(g(p,d)W ) > t => g(p, d) is separable. (20) 

Furthermore: The r defined in (|19J) is the minimal r with the property (|2C)|> . i.e. 

for all < t' < t there exists an entangled state g(p, d) with Tr(g(p, d)Wo) = t'. 
Proof. The idea of the proof is as described above: First, we compute the 
parameter such that for q < g_ we have KP(q) C SP(q). Via 1)16(1 we arrive 
at r. Finally we show that for q > q^ we have KP{q) <f_ SP(q). 

KP is just like BP a ball in the hypcrplanc P. The radius of BP is deter- 
mined by r 2 (BP) = (l-p) 2 d 2 -\\g(q,0) - g(p,0)\\ 2 = (1 -p) 2 d 2 - (3/4){q-p) 2 . 
Maximizing this over all p yields 

r 2 (^) = ^(l-<?) 2 - (21) 
The ball XP has the squared radius 

r 2 {XP) = ^-\{q) 2 . (22) 

If we choose q small enough, it follows from XP(q) C SP(q) that KP(q) C 
XP(q) C SP(q). So we can determine q- by the equation r 2 (KP(q-)) = 
r 2 (XP(q_)), apply (fT?)|) and arrive at r. 

Why can we not use a smaller r? The reason is that for any q e [0, 1/3] the 
ball XP{q) contains at least one state of dS. This means that from KP{q) % 
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XP(q) it follows that KP(q) % SP(q). To see the existence of a state G dS in 
XP(q) hrst note that if we take our first Bell state \ip+) = (l/\/2)(|01) + |10)) =: 
(1/V%)(|0'0') + |1'1')) and another Bell state = ^(|00) + 1 1 1> ) we have: 

((il - - |0 + )(0 + |))H5 =-J <0. (23) 

We can transform the bases |0'0') — > 1 00) , |1'1') — ► |11), in a continuous way 
and thereby construct a continuous map 

7 : [0,2] -> M; * 7 (*) (24) 

with 7(0) = an< i 7(2) = \4 >+ )(4> + \ an d for all t, 7 (i) is a maximally 

entangled state. There must exist a to (and without loosing generality we can 
set to = 1) with 

((il-|^ + )(^ + |)|(il-7(l)))H5 = 0. (25) 

This equation implies that 7(1) G 7 7 (q 1 = 0). The map 7 induces another con- 
tinuous map by 

7/:[0;l]^S; t » V (t) = ~1 + ± 7 (t). (26) 

Since 7 (t) is a maximally entangled state rj(t) is a point on SS. Furthermore all 
T)(t) have a distance of 1/VT2 to the state 1/4, i.e. 77(4) € 5(1/4, 1/VT2). Since 
7/(0) G P(l/3) and 7/(1) G P(0), for any g G [0, 1/3], contains one 77 and 
therefore XP(q) contains a state on dS. □ 
If one knows the value of p one can sharpen the statement of the last propo- 
sition: 

Proposition 3. Let p be fixed, g(p, d) 

*:=i-J— 
4 24p 

Then we have: 

Tr(g(p, d)W ) > =>• £>(p, d) is separable. (28) 

Furthermore: The $ defined in l|27|) is the minimal $ with the property 128|l i.e. 
for all < -d' < 1? there exists an entangled state g(p, d) with Tr(g(p, d)Wo) = 
Proof. The proof is essentially the same as the proof of Proposition 2. One just 
has to replace K by B p d and therefore KP(q) by BP(q). □ 

3.2 Error estimates 

Consider the case that one has measured Tr(Wog(p, d)) = a with a £ [0, r]. For 
such expectation value one cannot make a sure decision whether g is entangled 
or not. Nevertheless one can make a decision, if one accepts to make some error. 
Here we want to estimate the probability of making an error. First we give some 
analytical bounds on the error, then we perform some numerical simulations to 
rate the error. 

For an estimate of an error, we have to make one further assumption on 
g(p, d) : We have to assume some probability distribution. We always assume 
that for all p and d, g(p, d) is uniformly distributed in B Pi d- 



G B Pt d and 

3p {l-pfd 2 
8 2p 



(27) 
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3.2.1 Analytical estimates 



Our analytical estimation scheme relies on the fact that all states in the ball 
B(t/4, l/\/l2) are separable. Thus we know at least some separable states 
yielding expectation values a € [0, r]. If we assume that g(p,d) is separable, 
this gives us an upper bound on the probability of the error. 

For an illustration of this idea let us first assume that we know the fixed 
value of p and we have measured Tr(Wog(p, d)) = a. Via 1|16|) we can compute 
the q such that g(p 7 d) e V{q). We know then that g(p, d) G BP(q). If we now 
assume that g(p, d) is separable, the probability of guessing right is given by the 
ratio of the volumes 

_ vol(gP( g )) 
r ~ vol(BP(q)) ' [ yj 

and the probability of being wrong is bounded by 

< p vo\(BP(q)) -vol(X P(q)) 
ew ~ Ew = ^BPW) ' (30) 

These terms can be calculated with the standard formulas for balls in higher 
dimensions: vo\{BP{q)) = (tt 7 /7!)((1 - pfd 2 - 3/4(p - q) 2 ) 1 and vol {XP (q)) = 
U 7 /7!)(l/12-3/V) 7 . 

Since we do not know the value of p we have to maximize E w over all p, then 
we arrive at 

E_ = sup E w = l- V \ ^4 ■ (31) 

pe[o,i] (d{a + \)) 

This function is plotted together with numerical results in Fig. 1. 

Let us mention that with the same method as above one can also estimate 
the error for other scenarios, for instance if one generally assumes that for 
Tr(Wog(p, d)) > g is separable without looking at the precise expectation 
value; we will not discuss further this direction here, however. 



3.2.2 Numerical calculations 

One can investigate the properties of the witness also with numerical calcula- 
tions. For this purpose one can generate random matrices |16j and compare the 
result of the witness with the PPT criterion, which is a necessary and sufficient 
criterion of separability for 2 x 2 systems We have generated a sample of 
50000 random matrices (in Hilbert-Schmidt norm) in the ball -B(l/4, 1/^/12). 
By scaling we have retrieved then a set of 50000 random matrices g(p, d) for all 
p and d. With these matrices we have obtained our numerical estimates. 

First, we investigated the error estimate with E- from Eq. (|31J) . For this 
purpose we computed e w and then e_ as the supremum of e w over all p. This 
is shown in Fig. 1. 

One may ask how big the error is, if one generally concludes from Tr(Wog(j>, d)) > 
that g(p,d) € S. We have estimated the error (maximized over all p). The 
result is given in Fig. 2. 



4 Applications for 3-qubit and N x N systems 

First, we consider three-qubit systems. For these systems two types of tripartite 
entanglement are known: The GHZ-states and the W-states ^7]. Also families 
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Figure 1: Comparison between the numerical estimate (e_) and the analytical 
bound (E— ) as functions of the expectation value. The three curves are for 
three different values of d. 



of witnesses for these states are known ^H] and we show how to decompose 
them locally. Then we say something about a decomposition of a projector in 
N x N systems. 

4.1 Three qubits 

Three qubits can be entangled in different ways: They might be separable, biscp- 
arable of fully tripartite entangled |18j . The genuine threepartite entanglement 
consists of two classes: The GHZ-class and the W-class |17) . 

Witnesses for detection of GHZ-type states and W-type states have also 
been constructed in ^Hj. Here we want to show that these witnesses can be 
decomposed within our scheme. 

For the GHZ-class a witness is given by 

W GHZ = ^1-\GHZ)(GHZ\, (32) 

where \GHZ) is a pure state of the GHZ-class: \GHZ) = 1/V2(|000) + 1 1 1 1) ) . 
If q is a mixed state with Tt(qWchz) < the state g belongs to the GHZ-class. 
A decomposition of Wqhz can be achieved with similar calculations as above. 
The result is: 

W GH z = o (5 ■ 1 ® 1 ® 1 - 1 ® a z ® a z - a z ® 1 <Z> a z - a z ® a z ® 1- 
8 

— <T X <E> a X ® (J X + &X ® <%) <Ty + (7y ® <Ty g) U X + <Jy ® (J X ® tTy ) . 

(33) 

This witness can be measured with five collective measurement settings: Alice, 
Bob and Charly have to perform correlated measurements in the z-z-z, x-x-x, 
x-y-y,y-x-y and the y-y-x direction. 

For the W-states several witnesses are known. One example is the operator 

W W1 = ^1-\W){WI (34) 
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Figure 2: Probability of making an error (maximized over all p) when assuming 
Tr(Wog(p, d)) > g(p, rf) e 5 as a function of d. 



where |W) is now a pure state of the W-class: |W) = 1/V3(|100) + 1010) + 1001). 
This witness detects states belonging to the W-class and the GHZ-class, i.e. it's 
expectation value is positive on all biseparable and fully separable states. A 
decomposition is given by 

W W i = ^(13 • 1 ® 1 ® 1 - X ® 1 ® 1 - 1 ® ® 1 - 1 ® 1 ® + 

+a z ® a z ® 1 + er z ® 1 ® <j z + 1 ® a z ® a z + 3 • a z ® cr z ® <r z ) 
1 , 

- — (1 ® £T X ® &X + 1 ® Cj/ ® CTj/ + CTz ® <?x ® fx + &Z ® CTj/ ® + 

a x ® 1 ® cr x + cr y ® 1 ® a v + a x ® o-j, ® a x + a y ® a z ® ctj, + 

ffi ® cr^ ® 1 + (Ty ® cr y ® 1 + a x ® tTx ® (7j + cr y ® cr y ® a z ) . (35) 

Although this decomposition is a little bit longer, only seven correlated mea- 
surements are necessary. 

Another witness for W-class states is given by 

W W 2 = ^1-\GHZ)(GHZ\, (36) 

This witness can be measured locally with the same decomposition as l|33(l . It 
also can serve for a detection of states of the type (1 — p)l/8 + p\W)(W\ 7 as 
explained in |18j . 

4.2 N x N systems 

Now we want to generalize our results to higher dimensions. First we consider 
N x N systems, and the end of this section we make some remarks about N x M 
systems. 

A witness for an entangled state g with a non positive partial transpose can 
be constructed just like in the two-qubit case. First, one computes one eigenvec- 
tor corresponding to one negative eigenvalue of g Tn . The partially transposed 
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projector onto this vector is an entanglement witness. We want to decompose 
such witnesses for NPT states in this section. We only look at projectors, the 
partial transposition can be performed later. 

Our discussion proceeds as follows: After explaining our notation we con- 
struct a decomposition of a projector onto a state with Schmidt rank I using 
about 21 measurements. This decomposition is a generalization of the decom- 
position for the two qubit case. It is not clear whether this decomposition is 
optimal. Then we derive a lower bound for the number of measurements needed 
if the Schmidt rank I is maximal. We show that if I = N at least I + 1 measure- 
ments are necessary. 

We first explain some notational and technical details (the reader should 
consult \HJ\ for more explanations). We denote the real vector space of all 
Hermitean operators on Ha by HSa- In this space one can use the orthogonal 
basis {l,Gf,i = 1...N 2 - 1} where the Gf are the traceless generators of the 
SU{N), normalized to Tr{Gf) = 1. For N = 2 they are the Pauli matrices, for 
N = 3 the Gell-Mann matrices, etc. We can define Gq := 1 and can expand 
every projector (and any other element of HSa) in this basis: 

N 2 -l 

10) (01- E f&t (37) 

i=0 

where the entries of the Bloch vector /j are real, /o = 1/N, and from the fact 
that \4>)(4>\ is a pure state it follows that 

N 2 -l 

E # = l-i- (38) 



i=l 



We sometimes write (/o, Jn 2 -i) = (fo> f) = /■ It is easy to see that an 
operator described by g is a projector onto a vector orthogonal to \4>) if and 
only if g fulfills and 



N 2 -l 

<f,g>:= E fi9i = -jj- (39) 



One can also expand any projector (as every operator) on H.a ®Hb as 

N 2 -l 

iv>^i=E A ^® G ? ( 4 °) 

since the Gf cg> G^ form a product basis of the space 7iS = HSa ® HS b ■ 

Before we show our decomposition please note that is easy to decompose 
any operator A G TiS into iV 2 local measurements. One can always write 

N 2 -l N 2 -l /n 2 -1 \ 

A= E ^Gt ® Gf = E G?® E ^ G f (41) 

i,j=0 i=0 y j=Q J 

to obtain such a decomposition. This is also a decomposition of the form |J3}. 
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Theorem 1. Let \tf>) be a projector onto a state with Schmidt rank /. 
If I is even, \ip)(tp\ can be decomposed into 21 — 1 local measurements. If I is 
odd, can be decomposed into 21 local measurements. 

Proof. If we have = J2\=i s i|") we can write 

h/>}(^| = $>?R^| + E s i s Mh3) (42) 

i=l i,j=l,i<j 

with K(i,j) = + The first sum corresponds to one measurement 

and every of the 1(1 — l)/2 terms of the second sum can be decomposed by 
defining for every K(i,j) the directions \X^) = -±=(\i)±\j), \Y^) = -±=(\i)±i\j) 
and writing: 

K(i,j) = \X+X+)(X+X+\ + \X^Xr.)(X^X^\- 

-\ Y ijKj)( Y ij Y ij\ - \ Y u Y u)( Y u Y u\)> ( 43 ) 

as we have done before for 2x2 systems. This corresponds to 2 measurements 
for each K(i,j). 

The idea is now to sum up the terms from H43|l for different K(i,j) and 
K(m,n) in a way that the terms from different K(i,j) and K(m,n) can be 
measured with one measurement. 

Let us first consider the case that I is even. We have 1(1 — l)/2 index pairs 
(i,j). These pairs can be grouped into I — 1 sets of 1/2 pairs in a way that in 
every set every index 1 < i < N appears exactly in one pair. For instance for 
I = 4 the 3 sets may be defined as {(1, 2), (3, 4)}, {(1, 3), (2, 4)}, {(1, 4), (2, 3)}. If 
we look at the 1/2 K(i,j) belonging to one set, the corresponding vectors \Xf-) 
arc mutually orthogonal, they form an orthogonal basis of Ha and Hb- So all 
these vectors can be viewed as eigenvectors of some Hermitean operator on Ha 
and Hb and can be measured with one measurement. The vectors \Yfj) can 
also be measured with one measurement. So we need 2 measurements for one 
set and 2(1 — 1) measurements for all K(i,j). Finally we need one measurement 
for the first sum on the rhs of l|42() and this completes the proof for even /. 

If I is odd, we can similarly group the 1(1 — l)/2 index pairs into I sets of 
(I — l)/2 pairs. This time in every set every index appears at most one time, 
one index is missing in every set and every index is missing in exactly one set. 
As before we need 2 measurements for one set and therefore 21 measurements 
for all K(i,j). For the first sum on the rhs of (|42|l we do not need another 
measurement since we can put vector \i) to the set of index pairs where i is 
missing. □ 

Now we want to give a lower bound for the number of required measure- 
ments for a projector. This bound is based on the same idea as the proof of 
Proposition 2 and needs three lemmata. In Lemma 1 we give a lower bound for 
the rank of some matrix of the form (|14fl in N x TV-systems. In the Lemmata 2 
and 3 we show that the matrix coming from one measurement has a low rank. 
Together this proves our bound. 

Lemma 1. If \ip) G Ha ® Hb has the full Schmidt rank N then the ma- 
trix (Ay ) in igOll has the full rank N 2 . 

Proof. First, notice that the rank of (\j) is independent of the choice of the 
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basis Gf ®Gf. If one has another basis Hf <g> Hf with Gf = J2i a uHf and 
Gf = J2 r bjrHr the new matrix of coefficients is given by (AJ r ) = J2i j a Ti^ijbjr 
and since the matrices (an) and (bj r ) have full rank the matrix (A^) has the 
same rank as (Ay). 

Now we simply construct an orthonormal product basis of HS where (A£ .) is 
diagonal and the diagonal elements do not vanish. Starting from the Schmidt- 
decomposition \ip) = s i\w) WG define on Ha, as well as on Hb '■ 

Pk = \k)(k\, l<k<N (44) 

Qjk = -^(|i)<*| + |*)01). ±<3<k<N (45) 
Rjk = -j=(\j)(k\-\k)(j\), l<j<k<N. (46) 

These N 2 operators form an orthonormal basis of HSa (resp. HSb), denoted 
by Hf (resp. Hf), and if one computes 

N 

A; s = £ s a sp(a\Hf\p)(a\Hf\P) (47) 

a, (3=1 

one can directly verify that (X' rs ) is in the basis Hf ® Hf diagonal and has the 
full rank. □ 

Corollary 1. If \ip) G Ha <8 Hb has the Schmidt rank I then the matrix 
(Xij) in (gni has the rank I 2 . 

Proof. The proof is essentially the same as the proof of Lemma 1. We can view 
as a vector in a I x Z-system. □ 

Corollary 2. If |^>) £ Ha ®Hb has the Schmidt rank I then a decomposi- 
tion in the sense of requires I 2 Hermitean operators for every party. 
Proof. If one would need less, this would be a direct contradiction to Lemma 1 
and Corollary 1 . Please note that we have already computed this decomposition 
- see (HU. □ 

Lemma 2. Let V\, ...,v r € R n be some vectors obeying the equations 

<v i ,v j >=C^0 Vi^j. (48) 

v r should be uniquely defined by v\, ...,v r -i and the equations l|48|l while V r -\ 
should not be uniquely defined by v\, v r -2 and the equations (|48|l . Then we 
have 

v r e Lin(vi, ...,v r -i) (49) 

and 

dim(Lin(i/i, v r )) = r — 1 (50) 

where Lin({?i, ...,v r ) denotes the linear subspace spanned by v\, ...,v r . 
Proof. We can split v r in two parts: 

V r = V r \\ + V r± , (51) 
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where u r ii G Lin.(«i, v r -i) and v r j_A- Lin(vi, v r -i)- Since v r is unique, it 
follows that v r ± = (otherwise v r = v r ii — v r j_ would be a different solution) 
and the first part of the statement is proven. The equality in (f 501) comes from 
the fact that v r -i is not unique. □ 

Lemma 3. Let M be one LvNM in the sense of Q) expanded in the Gf <8 G^ 
basis: 

N N 2 -l 

M = ]T <% \Ai){Ai\ <g) |B 3 -) (Bj | = IMi Gf ® Gf (52) 

Then the (A^ 2 — 1) x (iV 2 — 1) submatrix in the right bottom corner of the 
TV 2 x A^ 2 matrix (/%•) (called (nij) re d — i^ij)i,j=i,...N 2 -i) has the rank N — 1. 
Proof. We can write any of the projectors and (Bj | as Bloch vectors 

Ai and -Bj (resp. Ai and Bj) with the help of l|37() . Then we have 

N 

Mred=Y, CijilifiBj). (53) 
»,J=1 

This is a (iV 2 -l) x (AT 2 -1) matrix, since every {Ai) T Bj is a (jV 2 -l) x (AT 2 -1) 
matrix. The range of this matrix is spanned by the vectors (Ai) T . The vectors 
(Ai) T correspond to the vectors \Ai), and they obey relations of the form (|3"5)) . 
Furthermore, An is uniquely determined by A\, An-i, since |^4jv) is uniquely 
determined by \Ai), |Ajy-i). Thus, we can apply our Lemma 2, and the rank 

Of {^ij) re d IS N- 1. □ 

Theorem 2. Let |^) € 7^ ®Hb have full Schmidt rank JV > 1. Then a local 
measurement of the projector requires at least A^ + 1 measurements. 

Proof. If we look at |'0)(V'l m the form l|4(J|l the matrix A,j has, according to 
Lemma 1, the full rank A~ 2 , the reduced matrix (\ij) re d — j=i,..jv a — l has 

a rank of at least A^ 2 — 2. 

Since the matrix (nij) re d corresponding to a single LvNM has, according to 
Lemma 3, the rank A'- 1 we need at least (N 2 - 2)/{N- 1) = A r +l-l/(A r -l) 
measurements. This proves the statement for A^ > 3. For N = 2 please recall 
that we have already computed that the submatrix {\ij) re d is of rank A^ 2 — 1 = 3, 
not N 2 — 2. This proves the claim for the case N = 2. □ 

The question remains, which of these results remain valid for A~ x M-systems 
with M > N. The answer is simple: All results remain valid. Since the maximal 
Schmidt rank in a JV x M-system is N, Theorem 1 can be proven in just the 
same way. Also the arguments which led to Theorem 2 can be applied. 



5 Bound entangled states 

In Hilbert spaces with dimensions higher than 2x3, there exist entangled states 
with positive partial transpose, the bound entangled states [301|^]. For this kind 
of states no general operational entanglement criterion is known and thus even 
complete knowledge of the density matrix may not suffice to decide whether a 
state is entangled or not. There exists, however, an important class of bound 
entangled states, the so-called "edge" states [22] , for which the optimal witness 
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operators can be constructed explicitly. In situations where an experiment is 
aimed at the generation of an edge state our method of local decomposition of 
a witness provides, therefore, a genuine experimental test. 

A state S is called an edge state iff it cannot be represented as S = qS' + 
(1 — q)(T s , where < q < 1, a s is a separable state and 5' is a state with a 
positive partial transpose. In other words, for all product vectors |e,/) and 
e > 0, 5 — e|e,/)(e,/| is not a bound entangled state anymore. This implies 
that the edge states lie on the boundary between the bound entangled states 
and the entangled states with non positive partial transpose. They violate the 
range criterion [201 in an extremal sense, i.e. S is an entangled edge state with 
a positive partial transpose iff for all product vectors |e, /) € R(6), |e, /*) ^ 
R(S Tb ), where R(S) denotes the range of S. 

The generic form of an entanglement witness for such a state S is |23| 

W = W-el, (54) 

where 

W = (P + Q Ta ) (55) 

e = inf (eJ\W\eJ), (56) 

and P and Q denote the projectors onto the kernel of 6 and 8 Ta , respectively. 

In the following we construct witnesses following this method for three kinds 
of bound entangled edge states from the literature, unextcndablc product basis 
(UPB) states introduced by Bennett et al. [21], "chessboard" states from Brufi 
and Peres [25] and the bound entangled states in 2 x 4 dimensions introduced 
by P. Horodecki |2(Jj . and decompose them locally optimizing the number of 
projectors onto product states or the number of settings as explained above. 
Note that for this particular construction it is from an experimentalist's point 
of view natural to decompose and optimize W rather than W , because the term 
el does not require any special setting. 



5.1 UPB states in 3 x 3 dimensions 

The states 

|Vo) = ^|0>(|0>-|1», hA 2 ) = ^|2)(|l)-|2», 
|V>0 = ^=(|0>-|1»|2>, |^ 3 ) = -L(|1)-|2))|0), 
l^4> = i(|0) + |l> + |2»(|0> + |l> + |2» (57) 

form a UPB [23, i.e. they are orthogonal to each other and there is no other 
product vector orthogonal to all of them. Therefore the state 



kl-X>>M)' (58) 



Qupb 

which is the projection on the space orthogonal to that spanned by the UPB, 
does not contain any product state in its range. Furthermore, it has a positive 
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partial transpose due to the orthonormality of the states \ipi). Therefore, g UF 
is an entangled edge state with a positive partial transpose. 
The projectors P and Q are related to each other by 



Pi = QT A =^2\1>i){1>i\, (59) 

i=0 

therefore we skip Q Ta and write the witness as 

4 

Wupb =531^X^1 -el- (60) 

i=0 

Five measurements are necessary to measure this witness, one for each of the five 
projectors, since the UPB is constructed in such a way that no two projectors 
can be evaluated in the same basis. The main problem of this construction is 
to find e. An analytical bound obtained by Terhal J2E] gives 

^1(6-V30)(2-^)^ 0QQ1297 (61) 

Numerical analysis leads however to the much bigger value e ~ 0.02842. 

In this case it is also interesting to optimize the number of local projec- 
tion measurements needed for the measurement of the total witness lf UPB . We 
extend the set {\ip%), i = 0, . . . , 3} with the vectors 

m = ^|0)(|0) + |l», fo&7> = -^|2>(|l) + |2», 

W = ^(|0) + |1»|2), |V>8) = ^(|1) + |2»|0), 

\H = |H) (62) 

to an orthonormal basis which can be used to decompose the identity. Alto- 
gether we are then left with a pseudo-mixture containing 10 projectors. Denot- 
ing B 1 ={\Q), |1>, |2)}, B 2 = {(|0) - |1))/V2, |2), (|0) + |l))/^>, B 3 = {(|1) - 
|2))/V2,|0),(|l) + |2))/V2},and J B4 = {(|0)-|l))/V2,(|0) + |l) + |2))/V3,(|0) + 
1 1) — 2 1 2) ) /2} , we easily see that measurement of W UPB decomposed in this form 
requires 6 correlated settings for Alice and Bob: B1B2, B^By, B1B3, B3B1, B4B4, 
and B\B\, therefore it is unfeasible to decompose the whole witness from the 
number of settings point of view as noted above. By subtracting in Eq. I|54|) 
some positive operator / instead of 1, one can reduce the number of projectors 
in the decomposition of W^vn to 9 - this gives an ONP, since the number of 
terms in any ONP must be larger than or equal to the rank of the witness, 
which is equal to 9. The idea is to form I as a convex sum of projectors onto 
|"0i)i=o....,4 and onto 4 other product vectors that are obviously not orthogonal 
to the 5 UPB states, but can be chosen such that the set of the 9 vectors forms 
a basis. If we choose as the additional vectors |i/>i)i=4 7 the decomposition 
contains 9 projectors in 5 settings. The bound for e has to be adapted as 

, . (eJ\W\eJ) 

€ = mt ~i — 7771 TT- (6 3 ) 

\ej) (e,f\I\e,f) 
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Numerical analysis leads to a value of e' ~ 0.0311. Note that when the bound 
entangled state is affected by white noise, namely p p = p ■ p VPB + (1 — p)t/9, 
the witness given above is still suitable for the detection of entanglement. For 
the witness in Eq. Ipjfi|l. Tr(W UPB p p ) < when p > (1 - 9e/5). 



5.2 Chessboard states in 3 x 3 dimensions 

The states introduced in are constructed from 4 entangled vectors, 



Pcb 

\Vi) 
\V 2 ) 
\Vs) 
\Va) 



?:=i 



(m,0,s;0,n,0; 0,0,0) 
(0, a, 0; 6, 0, c; 0, 0, 0) 
(n*,0,0;0,-m*,0;t,0,0) 
(0,6*,0;0,-a*,0;0,d,0), 



(64) 



(65) 



where * denotes complex conjugation and N = 1/ ■ By choosing the 
phases of the \Vi) and the basis vectors, 6 of the parameters can be made real. 
Without loss of generality, t and s can be assumed to be complex. In matrix 
form, p c b can then be written as 



N 



( 


m 2 + n 2 





ms* 











nt* 





\ 







a 2 + b 2 











ac 





bd 







sm 





M 2 





sn 

























a 2 + b 2 





be 





—ad 













ns* 





m 2 + n 2 





-mt* 













ac 





cb 





c 2 













tn 











— tm 





\t\ 2 













bd 





—da 











d 2 





\ 


























o / 



(66) 

Brufi and Peres suggested two methods of ensuring that p c ^ is bound entangled. 
We will employ the first one, i.e. we demand that p c b = p^ , which is fulfilled 
for t = ad/m and s = ac/n real and implies that P = Q. The kernel of p c b (and 
of p^if) is spanned by the (non- normalized) vectors 



l*i> 
|fc 2 ) 

1*3) 
1*4) 
1*5) 



|22) 
,m 

-,o, 

n 



a c 



(0. 



ac 



r,0; 



; 0,1,0; 0,0,0) 
be 



b 2 



0,1; 0,0,0) 



'- + b 2: 

(- — ,0,-; 0,0,0; 1,0,0) 
mn c 

bd ad 
(°- -7^-05^^,0,0; 0,1,0). 



b 2 '' 



(67) 
(68) 

(69) 

(70) 

(71) 



In order to decompose the witness W in terms of projectors onto product states 
we first examine whether there are more product vectors in the kernel of p c b, 
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so we try to solve 



(ai|0) + a 2 |l) + a 3 |2» <8> (&i|0) + b 2 \l) + b 3 \2)) = j^: 



(72) 



When writing down the equations one can see that the Xi can be substituted by 
products of one Oj and one Then it is possible to solve the set of equations 
which is then linear in the parameters bk- This in turn gives two equations for 
the dj . The first solution is given by 



and with 
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the other solutions can be written as 



|e,/) = ai(l,±72 ,1 ,7i 1 ' 1 



J ®6 2 f 
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™2 0,1 

72 



1,T- 
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ac7 2 ' 



(76) 



The parameters a\ and 6 2 can be used to normalize the vectors. We found 6 
product vectors in the kernels. Of those vectors, 5 will be linearly independent 
in general. Since they do not form an orthonormal set, we cannot construct 
P and Q from them. However, the witness can also be constructed by using 
instead of the projector onto the kernel of p c b (p^) an operator P (Q) which 
is strictly positive on the range of the kernel of p c h {p^ ) . This will only affect 
the value of e. Here P = Q can be constructed by summing the projectors 
onto five linearly independent product vectors from the kernel of p c b- Since the 
vectors are real, we have in addition that P = Q Ta . Hence in general the "pre- 
witness" W can be decomposed into 5 projectors onto product vectors requiring 
5 settings to measure Tr(W c bp), one for each of the projectors. 



5.3 Horodecki states in 2 x 4 dimensions 

The positive operators introduced in |2()j can be written in matrix form as 
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where b € [0, 1]. For b = 0, 1, the matrix pb is separable, and bound entangled 
for all other values of b. In the following we assume b ^ 0, 1. The kernel of pb 
is spanned by the entangled vectors 



= -^=(1,0, 0,0; 0,-1, 0,0) 
\k 2 ) = 4=(0,1, 0,0; 0,0, -1,0) 
= (0,0,l,0;y,0,0,-l), 



(78) 
(79) 
(80) 



where y = y(T— 6)/(l + b), so any vector in the kernel can be represented as 

\k} = (A,B,C,0;yC,-A,-B,-C) (81) 

where A,B,C,D are complex parameters. For \k) to be a product vector, it 
must be of the form 

|e, /) = (r, a) ® (A', B' , C, D') = (r(A',B', C , D')\ s(A' , B\ C',D% (82) 

where r, s, A' , B' ,C are complex parameters. It can be readily checked that 
there is no possibility to write |fc) in this form, therefore there is no product 
vector in the kernel of pi,. On the other hand, the kernel of p^ B is spanned by 
the entangled vectors 



|*i) = -/|(0,0, 1,0; 0,-1, 0,0) 



\/2 



V2 + y [ 



(0,0, 0,1; 0,0, -1,0) 
= (0, 1,0,0; 1,0,0,2/), 



(83) 
(84) 
(85) 



and does not contain any product vector, either. Therefore, the decomposition 
of the witness from Eq. (|54|) in projectors onto product vectors is a rather 
tedious task. On the other hand, we can write down the witness W as in Eq. 
(15411 and then decompose it as in Eq. @ and 1411) 



3 15 3 

w = ^2 Wi 3 ai ® t j = X! a% ® ^ 

i=0 j=0 



(86) 



i=0 



in a straightforward manner. Here di and n, for i > 0, are the generators of 
the SU(2) and SU(4), respectively, while a = l 2 /2 and t = l 4 /4. 
The matrices f turn out to be 
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(89) 



(90) 

but Alice and Bob 



where c = l/(2 + y 2 ). The number of correlated settings is 4 
need only 3 different settings each to measure the witness. 

Summarizing, we have shown that for UPB states and "chessboard" states in 
3x3 systems, one needs 5 LvNMs to detect the witness, whereas for Horodecki 
states in 2 x 4 systems, 4 correlated measurements are necessary. 



6 Conclusion 

In this paper we have studied the problem of detection of entangled states using 
entanglement witnesses and few local measurements. Two optimization scenar- 
ios were discussed: the one corresponding to the decomposition of the witness 
operator into the optimal number of projectors on product vectors (ONP) and 
the one corresponding to the decomposition of the witness into the optimal num- 
ber of settings of detecting devices (ONS). Several exact results and estimates 
have been obtained concerning optimal detection strategies for NPT states on 
2x2 and N x M systems, entangled states in 3 qubit systems, and bound en- 
tangled states in 3 x 3 and 2x4 systems. Despite numerous results and progress 
in understanding this problem, the general question of finding ONP and ONS 
for an arbitrary witness operator remains open. 

We wish to thank K. Zyczkowski, I. Cirac, S. Haroche, S. Huelga, B. Kraus, 
and H. Wcinfurtcr for discussions. This work has been supported by DFG 
(Graduiertcnkollcg 282 and Schwcrpunkt "Quantcninformationsverarbcitung" ), 
the ESF-Programme PESC, and the EU IST-Programmc EQUIP. 

While finishing this paper we became aware of a recent preprint by A. Pit- 
tcngcr and M. Rubin |27| . where it was shown how, if N is prime, a projector 
onto a pure maximally entangled state with full Schmidt rank can be measured 
with N + 1 LvNMs. Thus, our bound from Theorem 2 can be reached for this 
case. 
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